Parity Broken Chiral Spin Dynamics in Ba3NbFe3Si20i4 
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The spin wave excitations emerging from the chiral helically modulated 120° magnetic order in 
a langasite Ba3NbFesSi20i4 enantiopure crystal were investigated by unpolarized and polarized 
inelastic neutron scattering. A dynamical fingerprint of the chiral ground state is obtained, sin- 
gularized by (i) spectral weight asymmetries answerable to the structural chirality and (ii) a full 
chirality of the spin correlations observed over the whole energy spectrum. The intrinsic chiral 
nature of the spin waves elementary excitations is shown in absence of macroscopic time reversal 
symmetry breaking. 

PACS numbers: 75.25.+z,77.84.-s,75.10.Hk 



Ubiquitous in nature, chirality is what distinguishes 
a phenomenon from its materialization in a mirror [l[. 
Condensed matter physics and especially magnetism pro- 
vide a rich playground to investigate this phenomenon 
as it appears naturally to describe the arrangement of 
magnetic moments in non collinear magnets. The chiral- 
ity accounts there for the sense of rotation of the spins 
on moving along an oriented line. The latter can be a 
straight line around which the spins form an helix, or 
an oriented loop, for instance a triangle, at the summits 
of which the spins are orientated at 120° of each other 
(see Fig. [1]). For such coplanar spin configurations, it 
can be uniquely characterized on a bond connecting two 
consecutive spins Si and Sj by the parity breaking vec- 
tor product Xij — Si A Sj , defining the magnetic vector 
chirality. As far as domains of opposite chiralities are 
equipopulated, the chirality is usually averaged to zero in 
centrosymmetric magnets. However, neutron scattering 
experiments evidenced a pure macroscopic chiral phase 
in non-centrosymmetric compounds 0, 0] ? as f° r instance 
the Ba 3 NbFe 3 Si2 0i4 langasite [|. 

While these studies focused on the ground state prop- 
erties, the way the associated excited states inherit from 
the chiral properties remains an open issue. In this letter, 
we address this intriguing question, investigating the spin 
dynamics in a Ba3NbFe3Si2 0i4 single crystal by means 
of polarized neutron scattering experiments. Our study 
provides evidence for a fully chiral spin excitation spec- 
trum over the whole energy range. Chiral dynamics was 
already evidenced in a variety of magnetic materials [ful- 
fill , but we emphasize that these studies were carried out 
under a magnetic field, revealing the chirality through a 
macroscopically breaking of the time reversal symmetry. 
In the Ba3NbFe3Si20i4 case, the chiral dynamics is ob- 
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FIG. 1: from left to right: the two triangular chiralities ca = 
dzl, the two helicities ch = ±1, and the magnetic exchange 
paths Ji to J5 in Ba3NbFesSi20i4. 



served in zero field, thus ascribable solely to the space 
inversion symmetry breaking. To our knowledge, this 
has never been reported so far. 

The langasite Ba3NbFe3Si2 0i4 presents static chirality 
properties 0, [H, [Hj hereafter briefly recalled. It crystal- 
lizes in the acentric space group P321 which confers the 
compound a structural chirality denoted er = +(— )1 
for right (left)-handedness. The magnetic carriers are the 
Fe 3+ ions (5 = 5/2, L = 0), arranged on small triangle 
units (trimers) whose centers form triangular lattices per- 
pendicular to the trigonal c axis (see Fig. [1]). A magnetic 
order sets in below the Neel temperature T/v — 27 K with 
spins lying in the (a, b) plane and with the same 120° con- 
figuration on each trimer. This arrangement is helically 
modulated along the c axis with the period 1/r « 7 (see 
Fig. 4 of ref. 0). This magnetic order exhibits a single 
chirality of the helical order (named helicity), de- 
noted ch, and a single triangular chirality, denoted 
eA (see Fig. [I]). We note that Xij — S 2 sin(27rr)e^c for 
consecutive spins along c and XI ao Xij = (3v // 3/2)6' 2 eAC 
for the spins of an oriented triangle, with e# and 6a equal 
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FIG. 2: Inelastic neutron scattering intensities at 1.6 K in 
the (&*, c*) scattering plane on the time-of- flight spectrome- 
ter IN5. Displayed are intensity maps at constant transfer en- 
ergy (energy-cuts) and intensity maps along (0, k, — 1 + r) and 
(0,-1,^) lines in the reciprocal space (Q-cuts). In the bot- 
tom right sketch, the green area and arrow show the probed 
reciprocal vectors on the IN5 and IN20 spectrometers. 
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FIG. 3: (a-d) Magnetic scattering S(Q,lu) (red) and chiral 
contribution C(Q,cj) (blue) as a function of the transfer en- 
ergy at selected scattering wavevectors (0, l,-£), extracted at 
1.5 K by longitudinal neutron polarimetry on the triple-axis 
spectrometer IN20. The spectral weight is only significant for 
the branch emerging from the — r magnetic satellites, (e) Dis- 
persion of the chiral scattering C(Q, uj) obtained by gathering 
the measurements at all L 



to +1(— 1) for the right (left )-handed sense of rotation. 
According to 

mm, 

polarized inelastic neutron scat- 
tering can probe the dynamical spin chirality. This ex- 
tends the concept of static chirality to the dynamics by 
following the orientation changes during coherent mo- 
tions of spins in space and time. Its statistical aver- 
age is proportional to the antisymmetric off-diagonal dy- 
namical magnetic susceptibility (Tg|, [ItJ- We thus con- 
ducted experiments on the triple-axes spectrometer IN20 
at the ILL in its polarized neutron set up, with polariz- 
ing Heusler crystals as monochromator and analyzer. We 
used the CRYOPAD device, to obtain strict zero-field 
environment at the sample position, and to prepare in- 
coming and outgoing neutron polarization independently. 
Ef =14.7 meV was kept fixed yielding an energy resolu- 
tion of c± 1 meV and second order contaminations were 
removed by a PG(002) filter. The crystal, aligned with 
(0, k,£) as scattering plane, was cooled down to 1.5 K in 
an Orange cryostat. Different cross sections were mea- 
sured, a f^, cr^i and where the subscripts indicate the 



incoming/outgoing neutron polarization parallel t or an- 
tiparallel | to the scattering vector Q. Additional cross 
sections were measured at selected points in Q and en- 
ergy transfers Huj, allowing to separate the magnetic sig- 
nal from background, and to verify that phonon and in- 
coherent scattering are negligible in the investigated re- 
gion, could therefore be used as background above 
^1.5 meV. The magnetic dynamical structure factor was 
then obtained by S(Q,u) = ^ / (M y (Q,0)M y (-Q,t) + 
M z (Q,0)M z (-Q,t))e-^dt = (a n + a it )/2 - a n and 
the chiral magnetic dynamical scattering by C(Q,uj) = 
^ J(M y (Q, 0)M z (-Q, t) - Mz(Q, 0)M y (-Q, t)) e - <w *dt 
= (a^i — a\f)/2. Here, M ViZ (Q,t) are the Fourier trans- 
formed spin components at time t perpendicular to Q 
within (y) or perpendicular (z) to the scattering plane. 

To get first a global overview of the spin waves, an 
experiment using unpolarized neutrons was performed on 
the time-of-flight spectrometer IN5 at the ILL on a single 
crystal in rotation around the vertical zone axis a [18[- 
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FIG. 4: Dynamical structure factor calculated in the linear 
approximation from the minimal model (see text) to compare 
with the measurements of Fig. [2j 



The incident wavelength was fixed to 4 A, the chopper 
speed to 16000 rpm, yielding an elastic energy resolution 
~ 0.1 meV. Standard corrections and a background (high 
temperature scan) subtraction were performed. The data 
were then reduced with the Horace suite software [l9[ to 
obtain the excitation spectra as function of (Q,cj). 

Fig. [2] summarizes our inelastic unpolarized neutron 
scattering results obtained on IN5. Two spin wave 
branches emerging from the d=r magnetic satellites are 
identified [io|. They form delicate arches, with differ- 
ent maximum energies (lower branch clearly visible on 
IN5, maximum of the upper branch only observable on 
IN20). One of the branches is gapped with a minimum 
at around 0.4 meV whereas the other branch is found un- 
gapped down to the resolution (0.1 meV). A first unusual 
observation is the difference of intensities of the excita- 
tions emerging from the ±r satellites associated to a node 
of the reciprocal lattice. This effect strongly depends on 
the considered reciprocal lattice node and is clearly vis- 
ible for instance around (0,-1,-1) in Fig. [2] This is a 
signature of the structural chirality as detailed below. 

Fig. [3] gathers our inelastic polarized neutrons scatter- 
ing results of IN20. Energy-scans at constant Q were per- 
formed for different £ values along (0, 1,£) with £ ranging 



FIG. 5: Top: computed magnetic scattering S(Q,uj) for the 
two structural chiralities ct ± 1. Below each of them, com- 
puted chiral contribution C(Q,u) for the 2 associated mag- 
netic ground states (e#,eA = £h£t)- The red frame points 
out the calculation agreeing with the experiment in Fig. [3]e. 



from -r = — - to 1.7 (see sketch in Fig. [2]). The mag- 
netic scattering S(Q,u) (red lines in Figs. E^t-d) con- 
firms the magnetic origin of the two modes, although 
they are not as well separated as on IN5, due to the 
lower energy resolution. By extracting the chiral contri- 
bution C(Q,cj) (in blue), it is found that the lower mode 
is achiral (C(Q, u)=0) whereas the upper mode has fi- 
nite chirality (see Fig. [3]i and e). In addition, C(Q,uj) 
is positive and has a tendency to become negative for 
negative £ (compare Fig. [3^ and c at low energy). This 
change of sign is due to the fact that the neutron probes 
the moment component perpendicular to Q, C(Q, uj) cor- 
responding then to the projection of the Fourier trans- 
formed dynamical chirality onto Q. This result actually 
perfectly reflects a globally unchanged chirality of the 
upper branch all the way up to the maximum of the dis- 
persion. At large £ values, C(Q, uo) equals the amplitude 
of S{Q,uj) (see Fig. [3]i), pointing at a full chirality of 
the upper branch without chirality mixing. 

To shed light on this remarkable result, a character- 
ization of the ground state can be done from a careful 
analysis of the spin waves. We propose a model, based 
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on mean- field calculations, that accounts qualitatively for 
the magnetic order [4[ (see Fig. [1]). The in plane 120° 
spin arrangement is first stabilized by antiferromagnetic 
intra-trimer J\ and inter-trimer interactions within 
the (a, b) planes. Next, three inter-plane interactions J3 
to J5, connect each moment to the 3 moments of the up- 
per/lower trimer along the c axis. The acentric structure 
imposes different J3 and J5, oppositely twisted around c. 
The dominant one drives the helical modulation accord- 
ing to the structural chirality (J5 > Js,Ja for er = — 1 
and J3 > J^Jh for ct = +1). The two weakest inter- 
actions allow tuning the periodicity r of the helix. This 
model implies that the magnetic and structural chirali- 
ties are related (ca = ^t^h)' for the investigated crys- 
tal, ct = — 1 is observed (strong J 5 ), imposing opposite 
senses of rotation for the helicity and the triangular chi- 
rality. We are left with two different solutions (out of 
four) for the magnetic chiralities (en — +1, ca = — 1) 
and (ch = —1, ca = +1)- The observed ultimate se- 
lection of a single magnetic chirality was proposed [4j to 
originate from the antisymmetric Dzyaloshinskii-Moriya 
(DM) exchange interaction D.0i x Sj) with D the DM 
vector, allowed in absence of inversion symmetry center 
between spins Si and Sj. It suffices to consider the DM 
interaction inside the trimer with the same DM vector 
along the c axis for the three bonds [2l|. This favors 
planar spin components and its sign selects a triangular 
chirality 6a, and hence a helicity ch since 6a = ct^h- 

Using the standard Holstein-Primakov formalism in 
the linear approximation [22j, we computed S(Q,u) and 
C{Q,uj) at zero temperature. As shown by the com- 
parison of Fig. [2] and HJ a good agreement between ex- 
periment and calculation is achieved with the exchange 
parameters (in meV) Ji = 0.85 ± 0.1, J 2 = 0.24 ± 0.05, 
J 3 = 0.053 ± 0.03, J 4 = 0.017 and J 5 = 0.24 ± 0.05 
and a DM vector along c of ~ l%|Ji|. The latter was 
checked to produce the lower branch gap and to select 
a triangular chirality [23]. J3, J4, J5 were constrained to 
fulfill the r = 1/7 conditions and the set of best param- 
eters yields a Curie- Weiss temperature of 191 K, within 
less than 10% of the values obtained from susceptibility 
measurements 0, Hif . Other models have been proposed 
from spin wave measurements in Ba3NbFe3Si2 0i4, but 
they fail in reproducing all the features revealed in our 
measurements [HI, [26|. Interestingly, as shown by the 
calculations in Fig. [5j S(Q,lj) reflects the structural chi- 
rality: the asymmetric spectral weight of the branches 
emerging from the +r and — r satellites is inverted for 
ct = ±1. On the other hand, for a given et : S(Q,ou) 
does not depend on ea- The two corresponding magnetic 
ground states yield however C{Q,uo) with opposite sign 
(see Fig. [5]). S(Q,u) and C(Q,lj) therefore provide with 
a strong dynamical fingerprint of the ground state chiral- 
ity. Finally, the calculated chiral scattering for (ch = +1, 
6a = —1) reproduces very well the measurements of Fig. 



02 its sign, the absence of chirality of the lower branch 
and the full chirality of the upper branch. The calculated 
correlation functions show that the spins components in- 
volved are along a and b for the upper branch, and only 
along c for the lower branch, which yields a zero spin 
cross product and explains the absence of chirality. 

In summary, the spin waves excitation modes emerg- 
ing from the totally chiral magnetic order in a struc- 
turally enantiopure single crystal of Ba3NbFe3Si2 0i4 are 
extremely unusual. A spectral asymmetry is evidenced 
by inelastic unpolarized neutron scattering as a finger- 
print of the crystal chirality. In addition, the use of polar- 
ized neutrons and polarization analysis led to the discov- 
ery that the spin waves with in-plane spin correlations are 
fully chiral at all energies. Chiral dynamics was predicted 
from spin waves calculation for a Heisenberg triangu- 
lar antiferromagnet with unbalanced domains of opposite 
triangular chirality [17[ . This condition is obviously real- 
ized with a single chirality domain in Ba3NbFe3Si20i4. 
However, the absence of chirality mixing of the spin dy- 
namics is not only observed close to the magnetic satel- 
lites as calculated [17[ but up to the maximum of the 
dispersion branches («4.8 meV). It finally will be em- 
phasized that the crystal was in a strictly zero magnetic 
field (even the Earth's magnetic field is excluded from 
it) during the experiment on IN20. This establishes the 
observation of a dynamical chirality unbiased by macro- 
scopic time-reversal symmetry breaking, but solely asso- 
ciated with the space inversion symmetry breaking. 
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